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Abstract 

Within the context of flux tube models, heavy quark fragmentation takes 
place through the breaking of flux tubes with the production of a (di)quark- 
anti(di)quark pair. It is found that the (di)quark produced are more likely 
to be found in an L z = state. This naturally leads to an supression of the 
polarization distribution parameters W3/2 an d wi for (DijD^) and (A c i, A* x ) 
production respectively. The corresponding parameter w\ for (S c , E*) pro- 
duction, however, is not suppressed, in agreement with the CLEO results but 
not the DELPHI one. Implications on the measurements of Aq polarizations 
are discussed. 
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The production of excited heavy hadrons has been studied in Ref. 0. It was observed 
that, due to heavy quark symmetry and parity conservation of the strong interaction, the 
relative production probabilities of the different helicity states of the light degrees of freedom 
can only depend on the absolute magnitude of j, the helicity of the light degrees of freedom 
along the production axis, but not on its sign. As a result, the relative production probabil- 
ities of the (Di, D2) system is controlled by a single parameter W3/2, which is defined to be 
the probability that j has its maximal value |. A similarly defined parameter w\ controls the 
production of the (E c , E*) system, while another parameter A describes the likelihood of the 
production of a spin-1 diquark instead of a spin-0 one, which translates into the probability 
of producing a E c or E* instead of a A c . The framework has been subsequently extended to 
describe the excited A resonances (k c i,h* cl ) 0, with two more parameters (B,wi) defined 
in analogy to (A, Wi) for the (E c , E*) system. 

It must be emphasized that Ref. |J and @ are parametrizations rather than predictions of 
the fragmentation processes in the sense that they did not attempt to predict (or explain) the 
experimental values of the w's. By measuring the angular distribution of the decay products, 
ARGUS [[IH and CLEO [|5],|5| have measured u>3/2- They found w 3 / 2 to be small (best fit w 3 / 2 
is —0.30, which is in the unphysical region; restricting to the physical region < w 3 / 2 < 1 
gives W3/2 ~ 0), meaning that transverse polarization is preferred to longitudinal. There is 
no theoretical understanding of why w 3 / 2 is so small. In the case of excited B c production, 
one can calculate w 3 / 2 by perturbative QCD [0|§, and the result w 3 / 2 = 29/114 ~ 0.24 
is indeed on the small side. But it is dangerous to simply carry the result over to the 
(-Di,/^) where non-perturbative QCD effects are dominant. In the baryon sector, there 
is no measurement for w\ yet, while the w\ measurements by CLEO and DELPHI yielded 
inconsistent results. DELPHlQ obtained a small value for w± (best fit tfi is —0.36, which 



The DELPHI analysis use production from Z decays, not Ec . By heavy quark symmetry, 
however, w\ should be the same for both cases. 
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is also unphysical; restricting to the physical region < w\ < 1 again gives w\ ~ 0) flJlO 



while the CLEO result is consistent with an isotropic polarization (wi = 0.71 ± 0.13 while 
an isotropic distribution gives w± = |) fTT|j. 



3/ 

Intuitively the helicity of the light degrees of freedom can be viewed as the sum of two 
different contributions. One is the spin of the "brown muck" S, which is for A type 
baryons, | for mesons and 1 for E type baryons. Then this "brown muck" may orbit around 
the heavy quark with orbital angular momentum L, giving an additional contribution to the 
helicity. We will see that, if heavy quark fragmentation can be understood as breaking of 
color flux tubes as suggested by the Lund models [Q, the Artru-Mennessier model |13| and 



the UCLA model |14|, then the orbital angular momentum naturally prefers a transverse 



polarization, explaining the smallness of W3/2 and predicts a small w\\ w\, on the other hand, 
is not required to be small, i.e., the CLEO numbers are prefered. 

Due to the non-abelian nature of QCD, the color field is expected to be confined into 
tube-like regions (flux tubes) by the tri-gluon coupling. The flux tubes have constant ten- 
sion and ends at colored objects like quarks or diquarks. Such a picture is supported by 
Regge phenomenology, quarkonium spectroscopy, bag models and lattice QCD calculations. 
For concreteness, let's study the process Z° — > cc and the subsequent fragmentation and 
hadronization. Just after the Z° decay, both quarks are in general very off-shell and will 
fragment by the emission of hard gluons, which is governed by perturbative QCD. Even- 
tually such gluon bremsstrahlung will bring the off-shell energy down to Aqcd scale, and 
non-perturbative QCD will be important. This is when the flux tube model become a rea- 
sonable description of the dynamics. In a coordinate system in which the c quark travel 
along the positive z-axis in the cc center of mass frame (Z° rest frame if the momenta car- 
ried away by the hard gluons are negligible), there will be a flux tube lying along the z-axis 
joining the two quarks. The flux tube will be characterized by a constant linear energy 
density (tension) k ~ 0.2(GeV) 2 , which leads to a linear potential between the quarks. For 
a long flux tube, it will be energetically favorable to break the flux tube by the production 
of a quark-antiquark pair (or a diquark-antidiquark pair) to shorten the flux tube and hence 



reduce the energy stored in the flux tube. This is the QCD analog of pair creation in a 
strong electric field and is completely non-perturbative in nature. 

If we ignored the finite thickness of the flux tube, the pair-created quark- ant iquark pair 
will be produced right on the z-axis, where the flux tube is. Moreover, due to the tension of 
the flux tube, the antiquark will be linked by the flux tube to the c quark and move towards 
the c quark, i.e., along the positive z direction. In general there will also be transverse 
momenta, which has a gaussian distribution centered at zero, but let us ignore that for a 
moment. Then all colored objects (quarks and flux tubes) are on the z-axis, and the system 
has a rotational symmetry about the z-axis. As a result, the ^-component of the orbital 
angular momentum L z is conserved. Since the system starts out with vanishing L z (nothing 
is orbiting), the final hadron must have L z = 0. Notice that L 2 of the final hadron is not 
necessarily equal to zero, as the system does not have a spherical symmetry and hence L 2 
is not conserved. 

An alternative way of seeing the same result is to consider the wave function of the 
antiquark in the relative momentum space?]. Since the light antiquark is moving towards to 
c quark, it means that the relative momentum p of the light antiquark with respect to the 
c quark is in the positive z direction. In other words, the wave function of the antiquark 
in the momentum space is a wave packet peaked at some point po on the positive z axis. 
We do not know the exact shape of the wave function, but the rotational symmetry about 
the z axis mandates that the wave function can depend only on p z and ^Jp 2 + p 2 , but not 
the azimuthal angle 6. As a result, the expectation value of L z = idg vanishes. Notice 
that our argument holds even if the spread of the wave function is large with respect to \po\ 
and the wave function is non-vanishing even for points off the z-axis. As long as the wave 
function is azimuthally symmetric, L z has to vanish. Hence we see that L z in heavy meson 
production vanishes if transverse momenta are negligible. The polarization of the intrinsic 



2 An similar argument in the relative position space also holds analogously. 
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spin of the antiquark S, on the other hand, is not constrained in any way as long as it is 
cancelled by that of the quark produced at the same time. Since the flux tube models do 
not have a preference over the orientation of S, we naturally assume it to be isotropic. The 
total helicity of the light degrees of freedom, then, is the sum of of L, under the constraint 
L z = 0, and S, with an isotropic distribution. 

For a heavy meson with orbital angular momentum L 2 , the possible helicity states for 
the light degrees of freedom range over j = —L — ~, — L + — |, L + ~. By the 

conservation of parity, the probability of finding the light degrees of freedom with helicity j 
is the same as that with helicity — j. We will define Wj the probability of finding the light 
degrees of freedom with helicity either j or —j. The sum of all these probabilities equals 
to unity, i.e., J2Wj = 1. In particular, for L — 1, W 3 / 2 = w 3 / 2 defined in Ref. [Q, and 
W1/2 = 1 — u>3/2- The analysis above suggests that W1/2 = 1 and all other W's vanish, for 
all values of L. In other words, the light degrees of freedom will be in the lowest helicity 
state, and be as transversely polarized as allowed by quantum mechanics. In particular, for 
the [D X ,DD system, w 3/2 = W 3/2 = seen by ARGUS and CLEO [|H|. 

Our formalism can be extend to describe the production of heavy baryons as well. Differ- 
ent flux tube models have different descriptions of the production of the diquark-antidiquark 
pairs. In the simplest models, the diquark appears as a single entity, and the analysis above 
can be adopted in a straightforward manner. In some other models, like the "popcorn 
model" |Tj|, the two quarks in the diquark are produced in stages. The analysis above 
will be invalidated if the first quark-antiquark pair moves off the z-axis before the second 
pair is created, as such off-axis configuration will break the azimuthal symmetry. It turns 
out that, however, the first quark-antiquark pair will instead slide along the flux tube as 
"curtain quarks" but not wander off the flux tube. Hence, even in these models, all colored 
objects still lie on the z-axis, azimuthal symmetry is preserved, and hence the resultant 
heavy baryon will also have L z = 0. 

Since S = for a A type diquark, the light degrees of freedom of a A type heavy baryon 
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with orbital angular momentum L 2 can have helicity j = —L, . . . , L. Define as before Wj 
as the probability of finding the light degrees of freedom with helicity either j or —j. Then 
the analysis above suggests that Wo = 1 and all other W's vanish. In particular, for the 
(Ac^A^) system, W\ = wi as defined in Ref. [0, and Wo = 1 — w\. Our analysis then 
predicts w\ to be small. On the other hand, S = 1 for a £ type diquark, and the light 
degrees of freedom of a £ type heavy baryon with orbital angular momentum L 2 can have 
helicity j = —L — 1, — L, . . . , L, L + 1. Since S is supposed to have an isotropic distribution, 
S z is equally likely to be found in the +1, 0, or —1 states. This gives Wo = 3, W\ = | 
and all the other W 1 s vanish. For the (E c , £*) system, which is not orbitally excited, our 
analysis suggests an isotropic distribution, i.e., w± = ^ for w\ defined as in Ref. |IJ. This is 
in agreement which the CLEO result [JTTJ| but not the DELPHI one PJlO[|. 

Due to heavy quark symmetry, our analysis is obviously also applicable to b quark frag- 
mentation as well. It can also be easily generalized to other excited heavy hadrons. By 
the L z = rule, the light degrees of freedom of excited heavy mesons will have \j\ = |, 
excited A type baryons \j\ = 0, and for excited £ type baryons, \j\ can either be or 1, 
with probabilities | and | respectively. One exception is the P-wave £ type baryon, with 
the spin of light degrees of freedom si = L + S = . Then there will be only one helicity 
state and Wo = 1 trivially. In general, however, care must be taken to apply our analysis 
to very excited heavy hadrons, as the flux tube models may cease to be good descriptions 
with high excitation energies. 

It is interesting to see the implication of our analysis on the measurements of heavy 
quark polarization in Z° — * cc,bb processes. As discussed in Ref. |IJ, in general (except 
one special case) all polarization information of the heavy quark is lost in the meson sector. 
On the other hand, since the "brown muck" of Aq is spinless, the polarization of the heavy 
quark should be retained in the baryon sector. This effect, however, is modified by the 
presence of secondary Aq's, i.e., those produced in decays of excited heavy baryons. Since 
the polarizations of these secondary Aq's do not necessarily align with the initial heavy 
quark, the overall polarization is diluted. It has been shown || that, if one only includes 
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secondary Ag's from the (Eg, Eg) and (A c i, A* x ) doublets, the polarization is diluted by the 
factor 

l + |(l + 4 Wl ) + f(l + 4w 1 ) m 
1 + A + ' 1 j 

The first, second and third terms in both the numerator and the denominator correspond 
to Ag produced directly, from (Eg, Eg) decays and and from (A c i, A* x ) decays respectively. 
With W\ — |, W\ — and both A and B assuming the default Lund value 0.45 ||12|| , it is 
found that V = 0.65. For comparison, V = 0.58, 0.72 and 0.79 for W\ — Wi — 0, | and 
1 respectively. In the standard model, one find that the b and c quarks produced by Z° 
decay is partially polarized with Vb = —0.94 and V c = —0.67.0 Hence our analysis predicts 
V\ b = VbV = —0.61 and V\ c = V C V = —0.44. Note that this is two standard deviations 
away from both the ALEPH result |BJ V Ab = -0.23+J^(stat.)^o?(syst.) and the recent 
DELPHI preliminary result [[Tj]] V\ b = — 0.08^Q2g(stat.)^Q'jg(syst.). This should not be 
interpreted as the failure of our analysis, since the ALEPH and DELPHI central values 
—0.23 and —0.08 correspond to V = 0.24 and 0.08 respectively, which are not achievable 
for any choice of Wi and W\ anyway. Instead, this probably means that effects of other 
resonances, like the P-wave Eg's and D-wave Ag's, are not negligible. Also the hypothesis 
that A = B = 0.45 has not yet been tested. In fact, DELPHI preliminary results (1 < A < 2 
with large uncertainty) (H| suggests that more Eg'' are produced than the Lund model 
expects, and hence the depolarization is more severe. Evidently more accurate measurements 
on the parameters A, B and the u>'s are necessary to clarify the situation. 

The leading correction to our analysis comes from the transverse momenta acquired 
by the ( di ) quark- ant i(di) quark pair. These transverse momenta have a random (gaussian) 
distribution with (p±) = and (pj_) ~ (0.3GeV) 2 . Our analysis is a good description only 
if (p 2 z ) ^> (p±), i-e., the produced (di)quark-anti(di) quark pair move essentially in the z 



3 The definition of polarization used here differs from that in Ref. |l|] by a negative sign to facilitate 
comparison with the ALEPH results. 
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direction to preserve the azimuthal symmetry. Since both of these are governed by non- 
perturbative QCD, they should both be of order Aq C d, and it is not obvious that (p 2 z ) 
should be much larger than (p 2 j_). Physically, however, we do expect (p 2 ,) = \(p\) in the 
absence of the flux tube, and hence, under the tension of the flux tube, (pi) to be at least as 
large as |(pj_)- Consequently, we expect the probability of having orbitally excited hadrons 
with non-vanishing L z to be suppressed by 0(^(p 2 ± )/ (p 2 z ))- 

We conclude that, in flux tube models, orbitally excited heavy hadrons tend to have 
L z = 0. This gives a natural explanation of the small observed value of W3/2, and a small 
value is predicted for w\ but not w±. A large value of w\ (for an isotropic distribution w\ = |) 
would be fatal to our scheme, while a small w\ will mean that there are physics not captured 
by the flux tube model to make the diquark spin S anisotropic. It is expected that the next 
round of experiments will resolve the controversy on w\ and possibly measure Wi as well, 
putting our prediction to test. Our analysis is qualitative in nature; our ignorance of the (p 2 z ) 
prevents us from making quantitative predictions. On the other hand, it is important to note 
that our analysis is entirely non-perturbative in nature. This complements the perturbative 
calculations in Ref. |7|,|8| and suggests that the suppression of high helicity states is a genuine 
consequence of QCD, not that of the perturbative approximation. 
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